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Abstract. 

A semiclassical approximation for an evolving density operator, driven by a "closed" 
hamiltonian operator and "open" markovian Lindblad operators, is obtained. The 
theory is based on the chord function, i.e. the Fourier transform of the Wigner function. 
It reduces to an exact solution of the Lindblad master equation if the hamiltonian 
operator is a quadratic function and the Lindblad operators are linear functions of 
positions and momenta. 

Initially, the semiclassical formulae for the case of hermitian Lindblad operators 
are reinterpreted in terms of a (real) double phase space, generated by an appropriate 
classical double Hamiltonian. An extra "open" term is added to the double 
Hamiltonian by the non-hermitian part of the Lindblad operators in the general case of 
dissipative markovian evolution. The particular case of generic hamiltonian operators, 
but linear dissipative Lindblad operators, is studied in more detail. A Liouville-type 
equivariance still holds for the corresponding classical evolution in double phase, but 
the centre subspace, which supports the Wigner function, is compressed, along with 
expansion of its conjugate subspace, which supports the chord function. 

Decoherence narrows the relevant region of double phase space to the neighborhood 
of a caustic for both the Wigner function and the chord function. This difficulty is 
avoided by a propagator in a mixed representation, so that a further "small-chord" 
approximation leads to a simple generalization of the quadratic theory for evolving 
Wigner functions. 
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1. Introduction 

The Lindblad master equation describes the general evolution for markovian open 
systems under the weakest possible constraints [1] (see also e.g. [ll [3]). Given the 
internal Hamiltonian, H, and the Lindblad operators, Lk, which account for the action 
of the random environment, the evolution of the density operator may be reduced to 
the canonical form. 



Liouville-Von Neumann equation appropriate for unitary evolution. 

A typical example is based on the Jaynes-Cummings model, which describes the 
interaction of a two-level atom with a single mode of the optical field in a cavity. 
The statistically independent arrival of atoms leads to the damped harmonic oscillator 
equation for the photon field. 



where we identify the pair of Lindblad operators as proportional to the anihilation 
operator d = (q + ip)/ \/2 and the creation operator d'^ = (q — ip) / \pl for photons in the 
field mode (see e.g. [HE] and references therein). Further examples, e.g. laser models 
and heavy ions conditions, are reviewed in [6]. 

It should be mentioned that the Lindblad master equation is open to several 
criticisms. It has been explicitly shown, in the case of the dampened harmonic oscillator, 
that there is no universal master equation describing exactly the evolution of the reduced 
system, i.e. one that is fully independent of the initial state [7]. It may further be argued 
that other choices for the master equation can be more appropriate [8]. Nonetheless, 
the theory of semigroups is the most obvious generalization of the unitary evolution of 
isolated systems. Beyond its proven usefulness, this singles out the particular Markovian 
structure as a fundamental subject for investigation. 

The case where the Lindblad operators are all self-adjoint has deserved special 
attention. It is known that the corresponding Lindblad equation describes decoherence, 
or dephasing, as well as diffusion, but no dissipation [3] . Since the latter is usually a much 
slower process, it is often useful to simplify the evolution by considering only the self- 
adjoint part of the L^'s when studying the decoherence process (as in the semiclassical 
theory proposed in [9j). However, most physical processes for an open system such 
as ( 11. 2p are dissipative. It is therefore desirable to develop a semiclassical theory for 
the evolution of the density operator that combines the description of both the initial 
decoherence process and the more classical development of diffusion and dissipation. 

In this paper, we develop a formalism for treating the semiclassical limit of evolving 
density operators subject to equation fILlI) . including the cases where the Lindblad 
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operators are not self-adjoint. By semiclassical, we mean generalized WKB expansions 
(see e.g. [lO]), as opposed to simple power expansions in H. The present theory expands 
on our phase-space treatment for the semiclassical evolution of closed systems [Tl] and 
of non-dissipative open markovian systems [9] . To this purpose, we will adapt the theory 
developed in these papers (particularly [9]) in two respects. 

Recalling that R^^ stands for a (2A^)-dimensional phase space, which is a 
symplectic vector space, {x = (p,q)}, first we switch from the Weyl representation, 
where p is represented by the Wigner function W{x), to its Fourier transformed 
representation, the chord representation, where p is represented by the chord function, 
x(C); also known as the quantum characteristic function, given by 

^^^^ = (2^ / '^'^ ^^""^ exp { J(x A ^)} , (1.3) 
where we have used the skew product, 

N 

X A X' = ^ (Pn9n - InPn) = J X • x', (1.4) 

n=l 

which also defines the skew symplectic matrix J. The chord, £ = (^p,^^), is the Fourier 
conjugate variable of the centre x and stands for a tangent vector in phase space, as 
in the scheme for a Legendre transform. In contrast to the Wigner function, the chord 
function is not necessarily real, but its semiclassical expression is often similar to that 
of the Wigner function, as discussed in [I2l [T3] . 

In all cases where the Hamiltonian operator is at most quadratic in the momentum 
and position operators, x = (p, g), and the Lindblad operators are linear in x, as 
in example (11. 2p . the Lindblad equation reduces to a Fokker-Planck equation in the 
chord representation, which can be solved exactly for any initial state [H]. Various 
instances of this result have been previously reported, e.g. [IHl [lEl [TTj . Keeping to linear 
(but not self-adjoint) Lindblad operators, we now obtain an appropriate semiclassical 
generalization to the evolution of the chord function, x(^;^); for generic Hamiltonian 
operators, given an initial pure state, x(€)0)- This is similar to the theory for the 
evolution of the Wigner function in [9] , in which the Lindblad operators were assumed to 
be self-adjoint (no dissipation). But the present treatment has the immediate advantage 
of being exact in the quadratic case. 

In fact, the inverse Fourier transform of the semiclassical evolution for the chord 
function, evaluated within the stationary phase approximation, produces the same 
semiclassical evolution for the Wigner function as was presented in |[9j. However, this 
can now be seen to be a poorer approximation than the semiclassical chord function 
presented here. Indeed, the theory in [9] does not describe diffusion, which progressively 
coarse-grains the Wigner function. This is clear from the general analysis of the 
quadratic case [T^ , or [6l [18] in the case of initial coherent states. 

The second modification to the WKB semiclassical theory, which is required for 
treating markovian dissipation, is more profound: We find that it is necessary to work in 
double phase space, (x, ^) G R^^ x R^'^. This is a natural setting for the corresponding 
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description of the semiclassical evolution of the density operator [SU], or indeed, 
for the representation of general operators acting on the Hilbert space of quantum 
states, considered as superpositions of \ket){bra\ elements. Just as an evolving quantum 
state, \ip), corresponds to an evolving submanifold in simple phase space, x G R^^, the 
unitary evolution of a pure state density operator in a closed system, that is, a projector, 
corresponds to the evolution of a submanifold in double phase space. (In both 
cases, the respective submanifold satisfies an appropriate lagrangian property, to be 
specified). We thus obtain a formal generalization of the WKB framework, where an 
approximate oscillating solution of the Schrodinger equation is built from a classically 
evolving lagrangian submanifold [IDl EH |22] • 

However, the restriction to a closed system (and hence unitary quantum evolution) 
imposes a severe limitation on the allowed form of the corresponding classical double 
Hamiltonian [23]. The crucial point here is that an additional term in the double 
Hamiltonian arises naturally, as a consequence of the semiclassical approximation for the 
open terms in the master equation. This new term, which is responsible for dissipation, 
depends exclusively on the Lindblad operators and cancels in the special case where these 
are self-adjoint. Moreover, the resulting description of the full semiclassical evolution 
again coincides with the exact solution of the master equation in the quadratic case. 

Our use of classical double phase is limited to the semiclassical approximation. No 
attempt has here been made to define a generalized quantum mechanics that would 
correspond to classical double phase space. This could lead to a fully quantum path 
integral for markovian systems, as an alternative to the one developed by Strunz |24j . 
His approach relies on the position representation, with the obervables defined in the 
Weyl representation. Several of the ingredients in our theory already appear in Strunz's 
path integral, though its semiclassical limit is expressed in terms of complex orbits, 
whereas we deal only with real phase space propagation. 

This paper is divided in three parts. In the initial sections 2-5, we review 
basic material and reformulate the semiclassical theory for closed evolution |Tl] and 
nondissipative open markovian evolution [9], within the chord representation, so as to 
perfectly fit the exact quadratic results in |14j . 

In the second part, the ingredients in the basic result, equation ( 15. 2p . are 
reinterpreted within the double phase space scenario. This leads to the identification of 
the dissipative hamiltonian, in section 6, and the consequent semiclassical treatment of 
dissipative markovian dynamics, in section 7 (for linear Lindblad operators). 

It turns out that the classical region of double phase space, to which decoherence 
drives the evolution, projects singularly as a caustic onto the subspaces where either the 
Wigner function, or the chord function are defined. For this reason, in the final part of 
this paper, section 8, our semiclassical theory is adapted to the evolving centre-chord 
propagator [25] , which takes an initial density operator, expressed as a Wigner function, 
into a final chord function, thus avoiding caustics for a finite time. This leads to a small 
chord approximation for the evolution of the Wigner function itself. 
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2. Review of the semiclassical theory for density operators 

The chord representation of an operator A on the Hilbert space L^(R'^) is defined via 
the decomposition of ^ as a linear (continuous) superposition of translation operators, 

T| = exp{J(^Ax)}, (2.1) 

also known as displacement operators. Each of these corresponds classically to a uniform 
translation of phase space points Xq G R^'^ by the vector ^ G R^^, that is: Xq i-^ Xq + 
In this way, 

and the expansion coefficient, a function on R^^, is the chord symbol of the operator A: 
^(C)=tr(f„^i). (2.3) 
The Fourier transform of the translation operators defines the reflection operators, 



2^ R. 



exp{-(xA^)}f^, (2.4) 



such that each of these corresponds classically to a refiection of phase space R^^ through 
the point x, that is xq i-^ 2x — xq. The same operator A can then be decomposed into 
a linear superposition of refiection operators 



thus defining the centre symbol or Weyl symbol of operator A [26] , 

A(x) = 2^tr {R^ A). (2.6) 
It follows that the centre and chord symbols are always related by Fourier transform: 

(2^/'^'''^^''^"''^^^^''^^^^' ^^-^^ 

^(x) = I d^ i(^) exp { ^(^ A x)} . (2.8) 

In the case of the density operator, p, it is convenient to normalize its chord symbol, 
so that we define the chord function as 

whose Fourier transform is the Wigner function, 

Wi^) = I d^exp{'-i^ A^)} X{^), (2.10) 

or alternatively [27| 

iy(x) = ^-i^tr(^.p). (2.11) 
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The expectation value of any operator A, defined as 

(A) = trip A), (2.12) 
can then be written, according to (12. 5p . 

{A) = 2^1 A(x) tr (p R^) =Jd^ A(x) (x) , (2.13) 

which justifies the Wigner function being dubbed a "quasi-probabihty" , even though it 
can be negative. The normahzation condition reads 

1 = tr p = J dx W{x) = (27rfi)^x(0) • (2.14) 

The Weyl representation and its Fourier transform have a long history. References 
[26| [271 128| [29| ISOj [31] |32| [33] develop many of its aspects, with unavoidable variations 
in notation and interpretation. Our presentation is largely based on the review [34j- 

Standard quantum mechanical treatments requires us to choose between 
representations based on conjugate variables. This is just as true for the centre and chord 
symbols, related by (12.71) and (12. 8p . as for the more familiar position and momentum 
representations. However, the WKB semiclassical treatment links the x variable and 
the ^ variable through the stationary phase approximation. Indeed, starting from the 
integral expressions of (12.111) or (12.91) . this stationary phase method replaces respectively 
an integral over x, or ^, by its integrand, evaluated at one or several points Xj, or ^j. 
Because of the Fourier relation (12.101) between the pair of representations, each chord 
^ is then associated with a discrete set of "centres" x - this denomination will become 
clear in the following - while each "centre" x specifies a discrete set of chords 

This correspondence is geometrically clear in the case of a pure state, = \ip){i'\, 
classically associated with a (quantized) lagrangian sub manifold, C^, in the simple phase 
space X G R^^, that is, an A^-dimensional submanifold with the property that 



<j) p ■ dq 



0, (2.15) 

'7 

for any reducible circuit 7 lying in (see, e.g. [351 (SS] [37], for more on symplectic 
manifolds and their lagrangian submanifolds) . Then, for every point x, one can draw 
a discrete set of chords of the submanifold such that = x^ — x J and x 
is the midpoint of [x~,x+]. Reciprocally, every vector ^ coincides with a discrete 
set of chords for C^, with their midpoints at xj. These are the basic elements for 
the construction of a WKB semiclassical theory of density operators using this pair of 
conjugate representations, as was first noticed by Berry [32] . 

More explicitly, the construction of chords from centres, or vice versa, is realized 
as follows: To determine the set of centres that are conjugate to a given chord, ^, for 
each p, first translate the whole lagrangian submanifold, C, by the vector — ^, then pick 
the set {x"j} of all points of intersection between C and the translated submanifold 
The midpoint of each straight line, between x+j = x~j + ^ and x~j, defines 
Xj(^) = x.~j + ^/2, the centre associated to ^ [121 120]. To determine set of chords 
associated to each centre x, first reflect C through x and pick the set {x^j} of all points 
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of intersection between C and the reflected submanifold C^- Then, each reflected pair 
of intersections defines a chord associated to x [381 1201, i-e. ^j(x) = x+j — x~j. 

Given p and the corresponding C, the simplest semiclassical approximation for the 
chord function x(^) relates an amplitude aj(^) and a phase to each of the centres 

Xj(^) above, so that [I2l 

= E e^'^^(^)/^ = E XM) , (2.16) 

in such a way that 

x,(0=J^- (2.17) 



^.(x) = -J^ ■ (2.19) 



Similarly, the simplest WKB semiclassical approximation for the Wigner function 

iy(x) = E «i(x) e'^^W/^ = E , (2.18) 

i i 

relates an amplitude aj(x) and a phase 5'j(x) to each of the chords Cj(x) above, in such 
a way that 

.as;, 

9x 

The phases crj(^) (or ^^(x)) are also specified geometrically by half the action (or 
symplectic area) of a circuit taken along the original submanifold C and closed along 
the translated submanifold (or the reflected submanifold C^)- [§l The fact that the 
possible chords associated to a given centre always come in pairs (±^j) guarantees that 
the semiclassical Wigner function is real, as it should be. There is no such restriction 
for the chord function, unless the manifold itself has a special symmetry |12j . 

This simplest semiclassical approximation for the chord and Wigner functions is 
valid far from caustics, which arise for arguments of the chord function whose associated 
centres coalesce, or for arguments of the Wigner function whose associated chords 
coalesce, respectively. Hence, caustics are related to points of tangency between C and 
or between £ and £x, respectively [321 [381 [I2l US]. For the Wigner function, this 
occurs whenever x approaches C, in which case every pair of associated chords coalesce 
at a null chord (however, C is not the only region of centre caustics, generically) . The 
null chord caustic is more severe for the chord function, because in this case the entire 
manifolds £ and coincide. Thus, all the points in £ are associated centres to the 
null chord. 

The amplitude of each term in the above semiclassical approximation depends on 
variables that are constant along £. Defining the initial quantum state as an eigenstate 
of commuting quantum operators, the corresponding lagrangian surface, C (an A^- 
dimensional torus, if it is compact) will be defined by action variables X„(x) in 

§ Further Maslov corrections [10] should be included in the phase of the WKB semiclassical Wigner 
functions 32J and chord functions [12]. These are semiclassically small and do not alter the geometric 
relations ([2T9)) and ([2T7)) . 
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involution, i.e. all the Poisson brackets {X„,X„/} = 0. Let us now define the transported 
action variables, 

J±=J„(x±)=X„(x±^/2), (2.20) 

which may be resolved into either a function of x, for fixed ^, or vice versa. Then, 
generally, {Z^,X^,} ^ and it is found that the amplitudes are 

a(x) = I det{J+, J-,}|-i/2 ^ (2.21) 

within an overall normalization constant. This determinant can be reexpressed in terms 
of the Jacobian between the centre or chord variable and the 2N variables (X+,X~,) |38j : 

I det = I det{J+,X-,}| = I det ^^^\. (2.22) 

Clearly, the amplitudes, aj{^) (or aj(x)), depend on the degree of transversality of 
the intersection between C and (or £ and £x) and so they diverge at caustics 
|32l[38l[l2l[l3l. 

It should be noted that the equality between the amplitudes in both representations, 
equations fl2.2ip and (12.221) . holds for a specific pair of points (x^,x+) on the torus and 
hence for a specific centre-chord pair. In the centre representation, the Poisson brackets 
are considered as functions of x and we define x^(x). For the chord representation, 
these same endpoints are a function of ^ and so are the above Poisson brackets. The 
index, j, for the branch of the chord function (or the Wigner function) has been ommited 
from (12.211) . because a specific centre-chord pair (x, ^) will be a particular member of 
a set {(x, ^j(x))} for the Wigner function and, generically, a member of another set 
{(xj/(^), ^)} for the chord function. 



3. Review of the semiclassical limit for unitary evolution 

A theory for the semiclassical limit of unitary evolution, appropriate to density operators 
or unitary operators in closed systems, has been established in both Weyl and chord 
representations [HI [23l ESj. It is worthwhile to adapt the deduction of phase space 
propagators in [25] for the needs of the foregoing theory. The starting point is the 
product formula for any pair of operators, BA, in the chord representation: 

mm = J ho m - o e^^^"^^ (3.1) 

(see e.g. 01]). Here, when dealing with products of operators, we abuse the notation and 
use (i?A)(^) to denote the chord symbol C{$,) of the operator C = 13 A and, similarly, 
(BA)(x) stands for the Weyl symbol C(x) []]. 

The problem is that we will work with the chord representation of p, though the 
Hamiltonian should be specified in the Weyl representation. This latter is indeed a 
smooth function, H{x.), exactly classical, or at least close to it within the order of fi^, 

II Sometimes, the Weyl symbol of BA is denoted by the star product B-kA, when B is the Weyl symbol 
of B and A is the Weyl symbol of A. 
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{Hp - pHm = J §^ [H{^' + m - - m] Xi^') eTi-'<-^\ (3.6) 



whereas its Fourier transform, H{$,), is highly singular. By defining the translation of 
an operator as 

Arj := frjAf^r), (3.2) 
whose chord representation is given by 

^r7(0=e^'^"^i(^), (3.3) 

while its Weyl representation reads 

Ar;(x) = A(x + r7), (3.4) 

the phase factor in (13.11) can be incorporated into an integral involving both 
representations. Then, using (12. 4p . (13.31) and (13.41) . we rewrite (13. ip as 

{BAm = ^^-Ip J d^d^ A(x' - il2) B{i') et-'"(^-^\ (3.5) 

In this way, we obtain the chord representation of the commutator between H and 
the evolving density operator, p(t), as the mixed integral, 

d$!did ^ ^^^^ ^/oM .^^'^(^-^') 

{2-nhy 

Here, we emphasize, i^(x) is the Weyl representation of which is a smooth function, 
so that ( 13. 6p can be integrated in the stationary phase approximation. In the special 
case where -ff(x) is a polynomial, we can perform the integrals in (13.60 exactly. For a 
quadratic Hamiltonian we thus re-derive [H] 

(27rft)-^ {Hp - pHm = zn {h{^), x(0}, (3-7) 

which emulates the familiar result that the Wigner function evolves classically, when 
the Hamiltonian is quadratic |31j . 

For general Hamiltonians, we now insert the semiclassical approximation (12.161) 
for xi^yt) ill (13. 6p . Because of the linearity of the evolution equation for the density 
operator, it can be decomposed into branches Pj{t), each evolving separately, as 
represented by one of the semiclassical components, Xj(^) t) iii (I2.16p . Then (13. 6p can be 
integrated by stationary phase, to yield the lowest order semiclassical approximation: 

(27r/i)-^ {Hp-pH)sc{i) 

= E«.(€)(^(J^ + |) - ^(J^ - |)) e-^'^'" (3.8) 

= (^(x,(0 + - ff(x,(0 - ^/2)) x,(0 • (3.9) 
i 

Thus, by comparing with the unitary part of the master equation ( 11. ip . we find that the 
classical chord action aj{$,,t) evolves according to the Hamilton- Jacobi equation ^0^: 

similarly to the evolution for the centre action 
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It must be remembered that in fl3.10p . as well as in (13.61) through fl3.9l) . the function 
i7(x) is the Weyl representation of the quantum hamiltonian operator, which will be 
either identical, or semiclassically close to the classical hamiltonian function. In the 
general case where this hamiltonian function is nonlinear, the resulting evolution of the 
chord action is a consequence of the classical motion x^(^,t) of both the initial chord 
tips, x^(^,0) = Xj(^) ±^/2, whereas neither the chord, $, itself, nor the corresponding 
centres, Xj(^), will generally follow their respective hamiltonian phase space trajectories 
[TT] . By working directly on the double phase space, as discussed in section 6, a new 
Hamiltonian function can be defined on this doubled space to take account of the motion 
of both chord tips in a single trajectory [23] . 

The above approximation for the unitary evolution of the chord function does not 
include the evolution of the amplitudes, aj(^) in ( 12.16p . which can be obtained by 
including the next order in h in the theory. Alternatively, we note that, up to the 
leading order, the evolution can be portrayed as resulting from the full classical motion, 
i.e. all the trajectories generated by the hamiltonian, H{x), which transports the entire 
lagrangian submanifold, C{t), and its neighborhood. Thus, each pair of points, x^(^,t) 
on C{t) defines an evolving chord, 

|^.(t)=x+(^,t)-x7(|,t) (3.11) 

and an evolving centre, 

%itt) = (x+(tt)+xT(^,t))/2 . (3.12) 

One should note that, here, x^{^,t) denotes the hamiltonian trajectories of x^(^,0) = 
x^(^) = Xj(^) ± ^/2. Therefore, ^j{t) is generally different from the hamiltonian 
trajectory ^j(t) of the initial chord, ^j(O) = ^^(0) = ^, unless the hamiltonian is 
quadratic. Similarly, x.j{^,t) is generally different from the hamiltonian trajectory 

x,(tt) of x,(^,o) =1,(^,0) = x,(o m- 

By reconstructing the chord function according to the semiclassical prescription 
(I2.16P at each instant, the same phase evolution is obtained as from the Hamilton- 
Jacobi equation (I3.10p . but now the evolution of the amplitudes will also be included, 
as long as we also allow the action variables X„(x^) in (I2.2ip to evolve according to 
X„(x^,t) = X„(x^(t)), where ^^(t) is the hamiltonian trajectory of x''^(0). Again, 
it must be stressed that this semiclassical evolution of the chord function (or the 
Wigner function), resulting from global classical motion together with the geometric 
reconstruction of the representation at each instant, can only be identified with Liouville 
evolution (i.e. the evolution obtained from the hamiltonian trajectory of the argument 
of the the chord function, or the Wigner function) if the Hamiltonian is quadratic |llj . 

4. Chord representation of the open interaction term 

We now address various integral representations of the chord symbol for the open 
interaction term. The starting point is the product rule in the chord representation 
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(4.3) 



^^^^^^^^ = / ^§|p-^(^')^(^'')^(^''') '^(^-^'-^'-D exp [|:(^ A ^' - r A r 

where, again, we abuse the notation and write {ABC){$,) for the chord symbol of ABC. 
The exponent in the integrand is here one of the many different expressions for the 
symplectic area of the quadrilateral with sides: — ^. Incorporating the phase 

factor for translation into the chord representation, as in f l3.3p . leads to the compact 
expression, 

d$,'dC 

Even if the Lindblad operators L are not observables, as in the optical example (2), their 
Weyl representation are smooth functions on phase space, i^(x), whereas their chord 
representation, L{^), are quite singular. Therefore, we again need a mixed product rule, 
where a pair of operators, A and C are expressed in the Weyl representation: 

To obtain the desired expression for the nonunitary term of the Lindblad equation 
(II. ip . the order of the operators is permuted, which leads to sign changes for translated 
operators (13.21) . given by v4y^(x) = A(x + 77) in the Weyl representation, so that 

(ipi' - \i^U - \pmi) = I X(i') exp [ix' A (S - S') 

{L(x'+ |)L(x'- |)--i|L(x'+ |)L(x'+ |)- + L(x'- |)L(x'- |)-|}. (4.4) 

Note that L(x'+ ^/2)* is the Weyl symbol of the operator translated by ^/2, which 
is not equal to the adjoint of L^^^. 

The exact formula (14.41) is at a par with the representation of the commutator 
(13. 6p . It is interesting that, although (14.41) represents products of three operators, the 
dimension of the integral is the same as in (13.60 . Thus, including the presence of an 
internal Hamiltonian and a single Lindblad operator, the exact equation of motion for 
the chord function is given by 

XW'.O exp [l|x' A {-Wx'-|)-//(x' + |)] 

+ [l(x' + ^)L(x' - - -fL(x' + -)L(x' + + L(x' - ^)L(x' - ^ 

L V -r 2^ V 2^ 2^ ^ 2^ ^ 2^ ^ 2^ ^ 2 

If there are more Lindblad operators in the master equation (II. ip . then one must sum 
over these in the integrand on the right hand side of (14.50 . We have not included this 
obvious extension, so as not to confuse this sum with the further sum over semiclassical 
branches in the following formulae. 

Far from caustics, one can evaluate (14. 4p approximately, by stationary phase, if i^(x) 
is assumed to be a smooth function, by inserting the semiclassical approximation for each 
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separate branch of the chord function fl2.16p as in the previous section. The stationary 
phase condition singles out ^' = ^ and x' = Xj(^), one of the centres associated to a 
geometrical chord $, of the classical submanifold C. The full semiclassical approximation 
is simply 

(27rfi)-^(LpLt - h^Lp - \pPL)sc{^) = E {^(x,(€) + mHM^) - 

-^{|^(x,(0 + m\' + l^(x,(0 - e-^^^y\ (4.6) 

In terms of the chord tips, x^(^) = Xj(^) ± ^/2, the semiclassical approximation to the 
chord representation of the open interaction term can be rewritten as 

(27rfi)-^(LpLt - h^Lp - \pPL)sc{i) = 

E { - hn^tm - H^jmi' + ^ ^"^{^(x;(0)^(x-(0)*}} xa^), (4.7) 

where Xm denotes the imaginary part and Xi(^) is a branch of the semiclassical chord 
function given by (12.161) . In the case of a linear function, 

L(x) = / ■ X = r ■ X + i Z" ■ X , (4.8) 

as in the optical example 01. 2p . the semiclassical approximation for the open interaction 
term simplifies to 

(2vrfi)-^(LpLt - ]pLp - \pL'^L)sc{i) = 

E A A ^ - ^-[{l' ■ + (Z" • 1)2]) xA^). (4.9) 

On the other hand, (14. 4p can be integrated exactly, for a linear Lindblad operator 
(or even if it is a polynomial). Then (14. 4p becomes 

{2nn)-^{LpL^--L^Lp--pL^Lm = n{l' M") [{I' ■ i? + H" ■ if] x(0 , (4.10) 

in agreement with [14j. Compared with (14. 9p . we find the same second term on the right 
hand side. If the Lindblad operator is self-adjoint, i.e. I" = 0, this will be the only term. 
In this case, it is easier to develop a semiclassical theory for evolution of the density 
operator, which becomes exact in the case that the Hamiltonian is quadratic. This will 
be pursued in the following section. The first term was shown to describe dissipation in 
the exact quadratic theory ^1]. Though dissipation cannot be included in a standard 
semiclassical theory, we will show that it is naturally accommodated within the double 
phase space formalism that is developed in later sections. 

5. Decoherence without dissipation 

In this section, all Lindblad operators are restricted to be self-adjoint, so that 
Xm{L(x^(^))L*(x~(^))} = 0, simplifying equation (14. 7p . As pointed out in the 
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introduction, this means that the system may be considered to be conservative, albeit 
open to a random environment. 

If we further ignore the internal hamiltonian motion, or, more reasonably, restrict 
analysis of the decoherence process to its first stages, then we can consider the action 
(Jj{$,) to be constant in time, while the semiclassical amplitude evolves as 

a.itt) = «,(^,0)exp{ - |:5:|L,(x+(0) -^fe(x7(0)|^} . (5.1) 

k 

Generally, the above equation implies a fast shrinking of the chord function to a 
progressively narrower neighborhood of the origin. According to the discussion in |14j . 
this accounts for a fast loss of quantum correlations. However, for those chord tips, x^, 
that lie on a level curve (or level surface) of one of the real functions, Lfc(x), this term 
will not contribute to the loss of amplitude. The condition for a chord not to decay at 
all is that its tips should lie on the intersection of level surfaces for all the functions, 
Lfc(x). 

The effect of the internal Hamiltonian, H, can be included by considering the 
limit in a process where we switch it on and off, while alternatively connecting and 
disconnecting the Lindblad interaction (opening and closing the system) [9] . This defines 
a periodic markovian system in the limit of small periods, as in the periodization of 
hamiltonian systems in [39]. In the limit of short periods, there results a generalized 
Trotter ansatz [40j. Both the tips of the chord, x^(^), will evolve classically as 
x^($,,t) according to the Hamilton- Jacobi equation (13.101) for a time r/2, implying in 
the temporal evolution of a given chord, as $,j{$,,t) and for the centre, Xj(^), the 
motion xj(^,t), according to equations ( 13. lip and (13.120 . Then, at each opening of the 
system for a further period of r/2, the amplitude evolves according to (15. ip . Naturally, 
one must multiply both the open and the closed terms of the Lindblad equation by 
a factor of two, to make up for the reduced time in which either of them acts. In 
the limit as r — of an infinite number of closing and opening operations, we obtain 
the full semiclassical evolution of the chord function in a region free of caustics, as 
X{^,i) = Y.jXji^^t), with 

xACt)=X,%tt) exp[^D{x+(-t),x7(-t)}], (5.2) 

where xf{—t) is short for x^(^, —t) and X°j(^,t) denotes the semiclassical propagation 
for a time t of the j-branch of the chord function for the corresponding closed system 
(with all Lk = 0). The decay in amplitude for each branch of the chord function is 
determined by the decoherence functional over trajectory pairs, 

D{x+(t),x-(t)}:=5: fdt'\Lk{^+it'))-L,{^-{t'))\' , (5.3) 
k ^0 

where x^(0) = x(^) ± ^/2. Hence, it is the pair of backward trajectories ending at a 
given pair of chord tips on C{t) that determine the decrease in amplitude. The square 

% One should note that the proof of the Trotter theorem does not require fuU groups, appropriate for 
unitary evolution, but also encompasses semigroups, as in the present case. 
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root of the decoherence functional is a kind of time dependent measure of distance 
between any pair of points (x"^,x~), as pointed out by Strunz [21]. 

Concerning the derivation of the above semiclassical expression (15.21) . note that 
the chord function entering into the master equation is here the semiclassical chord 
function of a pure state (I2.16p . which is valid away from chord caustics. But far from 
the origin (a chord caustic), the damping factor in (15.21) is a non-oscillatory function. 
Hence, in a first approximation, it may be considered as a new factor of the semiclassical 
amplitude aj(^) in (15.21) . even though the exponent is divided by ft. Thus, when this 
modified expression for the chord function is inserted into the master equation, it is 
still the chord action function from (12.161) that defines the semiclassical evolution of the 
decaying chord function (15.21) . as long as all pertinent integrals are computed via the real 
stationary phase method. Accordingly, an improvement to (15. 2p could in principle be 
obtained by computing all pertinent integrals via the complex steepest descent method. 
This improvement is at present being investigated. 

The simplest case is where the Lindblad operators are all linear functions of position 
and momenta, Lfc(x) = • x. Then expression (15.21) simplifies, because 

L,(x+(|, -t')) - Lfe(x7(|, -t'))f = \lk ■ -t')f , (5.4) 

where the explicit dependence on ^ is emphasized in the r.h.s. Generally the evolution 
of each results from the hamiltonian flow of the tips x^(^), so that the evolution is 
j-dependent. However, if the internal Hamiltonian is a homogeneous quadratic, then 
the evolution of the chord is just given by |14| 



« = J f (5.6) 

and is therefore j-independent. Furthermore, the internal dynamics of the chord 
function is then Liouvillian: X^ji^^t) = Xj(^(~^))0)- this way, all j-branches of 
the semiclassical chord function can be combined into a single evolution, so that 

X(l, t) = x'i^i-t), 0) exp { - 1 ^ dt'lh ■ l(-t') I'}. (5.6) 

It is remarkable that this simple expression for the semiclassical evolution of 
an open system is actually exact and valid for any initial chord function (pure or 
mixed), under the above hypothesis for L and H [14J. Thus, no matter how full 
of quantum correlations the initial state might be, the infinite product of gaussian 
exponentials in (15.61) . or the more general exponential of the decoherence functional 
in (15.21) progressively squeezes them out. This process, by which the large chords are 
quenched, proceeds irreversibly, since D{x+(— t), x~(— t)} is a nondecreasing function of 
time. The semiclassical expression (15.21) generalizes the simple exact solution (15.61) . when 
the Hamiltonian is not quadratic, for chords that never lie close to caustics throughout 
the evolution. We retain the qualitative picture in which the evolving chord function is 
squeezed onto the origin by the decoherence functional, although D{x^(t), x7(t)} is no 
longer a quadratic function of ^. 
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The only possibility for the decoherence functional not to increase arises for pairs 
of classical trajectories generated by H{x.), lying along a level submanifold of the linear 
Lindblad-Weyl function, -L(x), i.e., the condition is that the Poisson bracket {L, H} = 0, 
which holds when the operators L and H commute. For more than one Lindblad 
operator, there is no dampening for those classical trajectories lying on the intersection 
of all Lfc level submanifolds, that is, when {Lfc,if} = 0, for all k. In the quadratic 
case, the specific evolution for each kind of classical dynamics (elliptic, parabolic or 
hyperbohc) is studied in [H]. 

The form of the evolution (15.21) goes some way towards justifying the rough 
qualitative description of decoherence (15. ip . that neglects the internal dynamics for 
very short times, since the chord function is seen to decay exponentially fast in the 
domain where ( 15. 2p is valid (which excludes a neighborhood of the origin). Generically, 
{Lk,H} 7^ 0, so that one needs (15. 2p to depict this quenching of the long chords in a 
fully quantitative manner. The further evolution of the state, for longer times, depends 
entirely on the remaining small chords, so that even our fuller semiclassical description 
is inappropriate. A generalization of the quadratic case (15.60 into the region of small 
chords is achieved indirectly in section 8. 

The semiclassical evolution of the Wigner function for open conservative markovian 
systems is obtained by the Fourier transform of the semiclassical evolution of the chord 
function. Each term of the sum, W^(x, t) = Y^jWj{'x,t), is given by a convolution 
integral of the unitarily evolving branch of the semiclassical Wigner function (12.180 with 
the Fourier transform of the decaying amplitude term in (15. 2p . This diffusive window, 
which coarse-grains the Wigner function, will broaden with time, as its inverse Fourier 
transform narrows down the range of the chord function. In the case of a quadratic 
Hamiltonian, the window will be gaussian and this description of the evolution of the 
Wigner function becomes exact [14J. 

The the evolving chord function (15.20 can only be inserted into the Fourier transform 
(I2.10p for chords that are far from caustics, which precludes small chords. For large 
chords, the convolution integral can be evaluated by (real) stationary phase, because 
the decaying amplitude term is a smooth function of ^ far from the origin. We then 
obtain a superposition of terms of the same form as ( 12.180 . each of them corresponding 
to a different branch of the centre action function Sj{'x,t). However, as with the chord 
function, the amplitude aj (x) now acquires a new time-dependent factor, so that we have 
the complete analogue of equation (15.20 for the semiclassical evolution of the Wigner 
function as 



where W^j{'x,t) denotes semiclassical propagation for a time t of the j-branch of the 
Wigner function as a closed system. The amplitude also decays according to the 
decoherence functional (15.31) that quenches the contribution of long chords. However, 
this is now determined by the choice of centre, x, rather than the chord, ^, i.e. here 
the chord tips are x^(x, t). One should note that the inclusion of the new decoherence 



iy,(x,t) = iy°,(x,t) exp D{x+(-t),xT(-t)} , 



(5.7) 
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damping factors in the semiclassical amplitudes of the chord function fl5.2p and the 
Wigner function (15. 7p does not destroy their equality (12.211) . 

The semiclassical expression (15. 7p was derived earlier in [9], solely within the Weyl 
representation. However, it does not lead to the exact evolution of the Wigner function, 
in the simple case of linear Lindblad and quadratic Hamiltonian operators, as does its 
chord similar (15. 2p . It follows that the theory for markovian open systems does not 
generalize the customary semiclassical covariance, among various representations of a 
given unitary evolution, with respect to Fourier transforms performed by stationary 
phase. Here, the chord representation has a decided advantage, because an essential 
qualitative feature is destroyed if the Fourier transform, leading to the Wigner function, 
is approximated by stationary phase. 

Thus, although ( 15. 7p well describes the initial stages of decoherence, it fails to 
address the diffusive process that sets in at longer time scales. In fact, there is no 
suggestion of the decoherence time threshold at which the initial pure-state Wigner 
function becomes positive definite, as always happens in the case of linear Lindblad and 
quadratic Hamiltonian operators [ITl [Ti] . The specific case of the damped harmonic 
oscillator is further analysed in [12]. This is the time which takes an initial state, 
represented by a Dirac-delta function in phase space, to evolve into a Gaussian with 
the width of a pure coherent state. At this time, any initial pure Wigner function 
evolves into a positive-definite phase space distribution, which is indistinguishable from 
a Husimi function [l3l HI] . In the general case of nonquadratic Hamiltonians, we can 
still define a local decoherence time as that which it takes the quenching factor in (15.20 
to shrink to the extent that it has the same area (or volume) as a coherent state. Beyond 
this time, its Fourier transform will coarse-grain away the fine oscillations of the Wigner 
function. 

6. Semiclassical evolution in double phase space 

WKB theory and its generalization to higher dimensions flUi [22] relates the solution, 
{q\ip{t)), of a Schrodinger equation to the corresponding classical Hamiltonian 
trajectories in the phase space x = {p,q) G R^^. This Schrodinger solution is the |g) 
representation of the unitarily evolving state, associated to a lagrangian submanifold 
C^, which is described within the Weyl formalism in section 3. This submanifold is 
more commonly described (locally) as a graph of the classical function, p(q) = ^{q), 
which maps q G onto p G R^ (which are also lagrangian coordinate subspaces, 
satisfying (I2.15P ). The classical action S{q) is both the generating function for and 
the oscillating phase of the quantum wave 

Analogously, the linear operators. A, that act on the quantum Hilbert space, form a 
vector space 1^4)), for which the dyadic operators \Q)) = \q~){q~^\ constitute a complete 
basis. Thus, defining the Hilbert- Schmidt product: 

tTA^B={{A\B)), (6.1) 
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(6.2) 



in close analogy to a wave function. □ It is then natural to relate a double Hilbert space of 
\ket){bra\ states to a double phase space: {X} = {x^} x {x+}, where = (p^, g^) (see 
e.g. [m 120] , or, for non- vectorial cases [35]). The operator \Q)) should then correspond 
to the lagrangian subspace, Q = constant, in the double phase space. This does hold, 
within a minor adaptation, due to the presence of the adjoint operator in the definition of 
the Hilbert-Schmidt product, or, more directly, the fact that, in ordinary Hilbert space, 
{bras\ are adjoint to \kets). Accordingly, if we define Q = {q~,q~^), we should define 
P = {—p~,p^) as conjugate coordinates on the double phase space {X = (P, Q)}. This 
is equivalent to changing the sign of the symplectic structure on R^^ = {x~}. 

In this way, we include, within the set of lagrangian submanifolds in double phase 
space, all the graphs of canonical transformations on single phase space, x^ i-^ x+ = 
C(x~). That is, we may rewrite the definition of a canonical transformation as 



where F is any curve defined on the (2A^)-dimensional submanifold, Ac, which is 
the graph of the canonical transformation C on the (2A^)-dimensional space {x~ = 
(g~,]9^)}, within the (4L)-dimensional double phase space, R^^ = {X = {P,Q)}. If 9 
is a parameter along F, then T{9) = (9) , •y^ (9)) , where 'y'{9) y^{9) = C{-y^{9)), 
and we may consider the curves 7^ as projections of the curve F. Going back to the 
operational meaning of this construction, if C~ is the lagrangian manifold corresponding 
to a quantum state | '?/'"), and C a canonical transformation, then = C{£~) can be 
interpreted as the lagrangian manifold of some 14''^) state, and the whole operation 
corresponds to a unitary quantum operator, Uc '■ ^ 

Besides portraying the graph of a canonical transformation as a Lagrangian 
submanifold, the product of a Lagrangian submanifold, C~ in {x^} with any another 
submanifold C'^ in {x+}, A = £^ x £+, is also Lagrangian in double phase space, 
but projects singularly onto either of the factor spaces {x"*^}. In the case that both 
submanifolds are tori, we obtain a double phase space torus, as if we had doubled the 
number of degrees of freedom. If = 1, it will be a 2-dimensional product torus [20] 
(taking care with the sign of p", in the present construction). 

If both Lagrangian submanifolds in single phase space correspond to the same state, 
i.e. = \ip~^), then we represent the corresponding pure state density operator, 

Pv — \'^){'^\ = l^))' the \Q)) representation as 



Therefore, its simplest semiclassical approximation can be expressed as a superposition 
of terms of the form 



It should be noted that we will not relie on the Hilbert-Schmidt norm and its evolution in the 
following discussion. 




(6.3) 




(6.4) 




(6.5) 
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with 

S,iQ) = Jq PAQ') ■ dQ' = fj ■ dg+ - J^'' pj ■ dq-. (6.6) 

Again, this is in strict analogy with the construction of semiclassical product states of 
higher degrees of freedom [2(1| . 

The next step is a change of lagrangian coordinates in double phase space: 

(P, Q) (x, y) , X = ^^-^ , y = J(x+ - X-) = J^. (6.7) 

Here, J is the constant symplectic matrix in single phase space and is used to canonize 
the initial tt/4 rotation on (x^,x+) that introduces the pair of lagrangian coordinates 
(x, ^) on double phase space. Thus, the pair of conjugate variables (x, y) also accounts 
for the sign change in the p~ coordinate. We should bear the discomfort that the 
canonical coordinate in double phase space is y, while the geometrically meaningful 
variable in single phase space is ^, the trajectory chord, which has x as its centre. It 
would also be possible to choose the variable, ^, as the conjugate to x, instead of y, but 
at the cost of writing the symplectic form on double phase space in a noncanonical way, 
leading to less familiar expressions for Hamilton's equations and some other elements of 
the semiclassical theory (see [15] for some of these expressions). 

If we consider the horizontal Lagrangian subspaces y = constant, each is identified 
with an element of the group of phase space translations, which includes the identity, 
since the identity subspace is defined as ^ = 0. On the other hand, the vertzca/ subspace, 
X = 0, defines the canonical reflection through the origin, x~ x+ = — x~ (or 
inversion), since all the chords for this transformation are centred on the origin (see 
[25] or [20] for further discussion.) 

We can now, in analogy to (16.61) . interpret the centre action 5'(x) in the semiclassical 
Wigner function (12.181) as 

5(x) = y(x') • dx' = ^ (x') A dx'. (6.8) 

The integral is evaluated along a path on the Lagrangian submanifold in double 
phase space, from some point on its intersection with the x-plane. (This intersection 
reproduces the single torus C^.) The integral is independent of the path on A^, because 
A^ is Lagrangian. We thus obtain the chord f l2.19p by taking the derivative of (16. 8p . 

The chord function is the Fourier transform of H^(x). If this transform of the 
semiclassical Wigner function is performed within the stationary phase approximation, 
the semiclassical expression for the chord function has a phase, cr{^)/ti, such that the 
chord action, c'"(^), is the Legendre transform of the centre action, 5'(x). It can be 
defined directly in terms of a similar integral to (16.81) . with the roles of x and ^ reversed: 

a(C) = x(^') A d^' = - f'^ x(y') ■ dy' = a'(y). (6.9) 

JO Jo 

The action cr(^) is, of course, the same as appeared in the semiclassical theory for 
the chord function (12.161) . When this theory is transported into double phase space. 
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it is often simpler to deal with o"'(y). Then, within this formalism, the semiclassical 
expression, for (each branch of) the pure state Wigner function or chord function, 
assumes a generalized WKB form, derived by Van Vleck [16] . 

So as to treat the unitary evolution of the density operator, which preserves the 
purity of the state, \'ip){ip\, we need to consider the corresponding classical evolution of 
both the tips of each chord, x~ and x"*", lying on a 2A^-dimensional lagrangian torus. 
Taking account of the sign change of p~ , in the definition of double phase space, we find 
that the double phase space Hamiltonian must be [23j 

MuiX) = if(x+) - H{^-) = - Jy/2) - /J(x + Jy/2). (6.10) 

This Hamiltonian dynamics evolves lagrangian submanifolds in double phase space, 
which correspond to pure density operators satisfying the Liouville-Von Neumann 
equation. The explicit formulae for the semiclassical evolution of the Wigner function 
are given in pT| l23]. whereas the evolving action (13.101) of the chord function is presented 
in [25]. Reinterpreted as an evolving action in double phase space, (IS.lOp assumes the 
form of an ordinary Hamilton- Jacobi equation for cr'(y): 

The difficulty lies in the neighbourhood of caustics of the initial state, which require 
more sophisticated semiclassical treatment. 

According to the discussion in section 3, the evolution of the amplitudes, in the 
decomposition of either the Wigner function or the chord function, relies on the previous 
specification of translated action variables, for the lagrangian manifold, C, corresponding 
to a semiclassical state. The corresponding lagrangian submanifold is now a product, 
A = X £+. Thus, the quantized double torus is defined as the intersection of all the 
level submanifolds of the 2N variables X^, defined by (12.201) . or 

T±(X)=X„(xTJy). (6.12) 

The fact that we are dealing with projection operators restricts the Bohr-level for each 
pair of variables (X;J"(X),X~(X)) to be the same. 

Let us then consider the family of actions, S'(x, X^) or a{^,I^), evolving classically 
for all possible constant values of the action variables, X^. This is known as the complete 
solution of the Hamilton- Jacobi equation [36j. Then, a simple extension to double phase 
space of the usual canonical formalism implies that 

,^ ^ gV(y,X^t) , 5X± 1 ,^ ^ g^g(x,X^,t) , aX± 1 

Idet^^^^| = |det^| and I det ^^^^^ I = I det ^ | . (6.13) 

Combining (16.131) with the expressions (12.211) and (12.221) for the semiclassical amplitudes, 
leads to 

, , a2^(x,X±,t),i/2 '9V(y,2:±,t),i/2 ^ 

= |det Q^QjJ ' \'^' = \det QyQ^^^ l'^' = «'(y,i) = ^i^t), (6.14) 

* If this condition is relaxed, the present semiclassical theory is immediately extended to include 
the propagation of dyadic operators, corresponding to pairs of different eigenstates of N commuting 
operators. 
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where this equahty fl6.14l) between the centre and the chord amphtudes holds only at 
a specific double phase space point X = (x, y) on the double torus A, as pointed out 
in the discussion at the end of section 2. Therefore, the amplitudes of the evolving 
Wigner and chord functions are entirely determined by the complete solution of the 
respective Hamilton- Jacobi equations, in full analogy to the ordinary evolution of WKB 
semiclassical states in the Schrodinger formalism, derived by Van Vleck [16]. 

It should be pointed out that our use of double phase concerns only the semiclassical 
approximation to the evolution generated by the Liouville-Von Neumann equation and 
susequently the full Lindblad equation. At each instant, either the Wigner function, 
or the chord function are defined in the standard way, as the traces (12. 6p . or (12. 9p . 
in terms of the single phase space of centres, x, or chords ^, respectively. Though it 
may be tempting to define an enlarged quantum evolution for superoperators in direct 
correspondence with double phase space, no such generalization is treated here. 

Now, we finally turn to the semiclassical theory for markovian evolution, as 
discussed in the previous section. We immediately recognize in (15. 3p the same structure 
as that of the double phase space Hamiltonian (I6.10p . that is: 



This double phase space Lindblad function is the basic ingredient in the decoherence 
functional, which is now defined along a single trajectory in double phase space, 
generated by IHu{X.): 



The square root of this functional can now be interpreted as a time dependent length of 
the double phase space vector, X(t), with D{(x, y = 0)} = for all time, instead of a 
distance between a pair of of single phase space points. 

In conclusion, we can interpret the conservative semiclassical evolution of the chord 
function (15. 2p entirely within the double phase space picture. Indeed, this has assumed 
the same form as general Van Vleck evolution, with the centre variables x playing the 
role of positions, while y = stand for the momenta. Following this analogy, the 
Wigner function substitutes the Schrodinger wave function and the chord function is its 
Fourier transform. 

The only new element that has been added is the action of the decoherence 
functional: The amplitude of the chord function away from the origin progressively 
decays in time. This quenching of the long chords can only be partially incorporated 
in the semiclassical Wigner function, obtained from the alternative projection of the 
same Lagrangian submanifold, A(t), in the limit of very short times. Indeed, it is only 
through this submanifold that we can ascribe specific chords to each centre x. For finite 
times, it is better to calculate the Wigner function as a full convolution, according to 
the discussion in the previous section. 

So far, that is, in the semiclassical theory for quantum unitary evolution, the 
concept of double phase space may be considered to be somewhat redundant, because 



iL(X) = L(x+) - L(x-) = L(x - Jy/2) - L(x + Jy/2). 



(6.15) 




k 



(6.16) 
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everything can be described in terms of pairs of hamiltonian trajectories in single phase 
space. However, and this is the fundamental point, for dissipative markovian systems, 
we can identify the dissipative term in the full semiclassical master equation (14. 7p . that 
is, 

1Hl{X) := - 5] Xm Lfc(x+)Lfe(x7)* = ^ Jm Lfc(x - Jy/2)Lfc(x + Jy/2)* , (6.17) 

k k 

as a new term of the double phase space Hamiltonian. Indeed, the introduction of (13. 9p . 
together with the generalized version of (14.71) for several Lindblad operators, into the 
chord representation of the master equation (11.11) results in the semiclassical evolution 
equation. 




-^E|^^(4)-^^(^7)| }xA^^t), (6.18) 

for each branch of the chord function, recalling that x^ = x^(^) = x ± ^ = x =l= Jy. 

In this way, for open dissipative systems, we can consider the total Hamiltonian 
function on double phase space R^^ x R^^ as given by 

1H{X) = MuiX) + MLiX) , (6.19) 

where MuCX.) and IHi{X.) are given respectively by equations (16.101) and (16.171) . for 
functions H{x.) and Lk{:x) on simple phase space R^^. The particular combination of 
Lindblad functions, which we have recognized in (I6.17P as a new term in the double 
Hamiltonian, also appears as the integrand of the phase functional in Strunz's path 
inegral [21], but there, in the absense of double phase space, it lacks an interpretation. 
Indeed, it is the very fact that 1Hl(X.) cannot be related to a Hamiltonian in simple 
phase space in the same way as ]Hu(X.), which now establishes the double phase space 
formalism as wholly indispensable. Therefore, for open dissipative systems, a trajectory 
of the full Hamiltonian 1H(K) in double phase space is not equivalent to a pair of 
trajectories of a Hamiltonian in simple phase space, as in the semiclassical theory for 
closed systems. 



7. Dissipative semiclassical evolution 

To simplify our study, we here restrict the Lindblad operators to be linear functions 
of positions and momenta, that is, given by (14.81) . The dissipative term (16.171) in the 
double phase space Hamiltonian is then rewritten as 

i?fi(X) = -7x-y = -7^Ax, (7.1) 

defining the dissipation coefficient, 

7 = E A II (7.2) 
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The contribution of this term to Hamilton's equations in double phase space is 

x=-7x , y = 7y , (^ = 7^), (7.3) 

so that the dissipative evolution in double phase space is always hyperbolic: contraction 
on the lagrangian subspace R^^ = {x} (the identity plane, i.e. the centre phase space) 
and expansion on its conjugate R^^ = {y} (and hence of the chord space), or vice versa. 
The fact that these are precisely the subspaces which support the Weyl and the chord 
representation singles them out as privileged choices for the description of quantum 
markovian processes. 

From here on, we unleash the centre space from its interpretation as single phase 
space. This aspect has been previously addressed in the review of Balazs and Jennings 
[33], but now it becomes an intrinsic feature of the present theory. For, it is by 
recognizing, instead, the centre space as a lagrangian subspace of the double phase space, 
that we free symplectic areas and volumes on the centre space from being conserved. In 
this way, dissipation finds a place in a strictly real Hamiltonian theory. |J 

For the case where L is the annihilation operator L = d, then /' = (0, 2~^/^) and 
/" = (2~^/^, 0), so 7 > and the centre motion is contractive, while the spacing between 
neighboring chords expands with time. In the case where L is the creation operator a^, 
the opposite happens. The optical master equation (11.21) combines both the creation and 
annihilation operators, but in such a way that 7 > 0, so the centre motion is contractive 
(strictly dissipative). In the example of a two level atom coupled to a bath of photons 
in a single field mode, this indicates that, although stimulated emission compensates 
absorption, spontaneous emission leads to an irreversible loss. 

The double phase space formulation maintains the correspondence of the evolving 
density operator, p(t), to a time-dependent lagrangian submanifold, A(t), even in the 
presence of dissipation. It is true that, unlike lHu{X.), the new dissipative term, 1Hl(X.), 
of the double phase space hamiltonian destroys the factorization of the initial double 
lagrangian submanifold into single phase spaces. However, the generalized Van Vleck 
form of the evolution, which was shown to hold for phases and amplitudes in the previous 
section, is in no way restricted to product tori. Thus, there is no obstacle to the 
immediate generalization of the present theory. 

The evolving lagrangian submanifold, A(t), can in principle be described by an 
action function that measures its symplectic area with respect to any (double phase 
space) lagrangian coordinate subspace, such as the position subspace, Q = {q'',q~^), 
even though we can no longer describe A(t) locally as = ± ds{q,t)/dq^, with 
the same action function, s{q), for q~ and q^. Nonetheless, the Weyl and the 
chord representations, in terms of the double phase space variables, x and $, = — Jy 
respectively, have privileged roles: The former continues to be interpreted as a 

However, the symplectic structure on the centre space remains indispensable, as concerns its role in 
the definition of the Weyl symbols. The evolution of this representation follows from the motion of the 
density operator, not from an evolution of the symplectic structure on the centre space, which is not 
preserved by the hamiltonian flow of (j6.10p . even if L = 0, as long as H is not quadratic. 
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quasiprobability, whereas the evolution of the latter exhibits decoherence and diffusion 
in a specially simple form, by merely quenching the amplitude of the various branches 
of the chord action function (16. 9p . 

The evolution of each branch is then obtained from the same Hamilton- Jacobi 
equation (16. lip as before, except that now we replace the closed double Hamiltonian 
(I6.10p by the full markovian double Hamiltonian (16.190 . Spelled out in terms of the 
single space functions, for 1Hl{X) given by (17. ip . this becomes 

dcr ,^ . ^^,^da £, ^^,^da £, ^ da ,„ 

= ^(J;^ + 1' + -a? • f^-*) 

Recalling the simple form of the action ( 16. 9p in double phase space, the evolution of the 
lagrangian submanifold, A(t), is just given by x(y,t) = — (9cr'(y, t)/9y. 

The semiclassical approximation for the evolving chord function is still given by 
(15. 2p for each of the branches of the chord function. The action function evolves 
according to (I7.4p . whereas the amplitude of the decoherenceless factor, ji.^^'t) is 
specified by (I6.14p . In the general case, where the Hamiltonian is not quadratic, 
the closed evolution of the decoherenceless factor of the chord function in (15. 2p is not 
obtained from the single phase space trajectories, ^(— t), generated by the classical 
Hamiltonian, if(x), i.e. ^ 0). 

Having reinterpreted the decoherence functional (15. Sp in double phase space as 
(I6.16p . we now obtain the full semiclassical markovian evolution for each branch of the 
chord function, including dissipation, in the same form (15.20 as before, even though the 
evolving ingredients can no longer be interpretated in single phase space. Because of 
the linearity assumed for all the Lindblad operators (14. 8p . the decoherence functional 
(I6.16P or (15.30 takes the explicit form, 

D{x;(t),xT(t)} = ^ f dt'[\l',.l^{i,t')\' + (7.5) 

for each branch of the chord function. It is the classical evolution of the full double 
phase space vector, Xj = (x~(^), x+(^)) (or, in alternative coordinates, = 
(xj(^), y = J^)) that determines the decoherence functional, i.e. Cj(C, —t') is obtained 
by multiplying the y-component of Xj(t) by — J. 

It is only the further restriction to a quadratic Hamiltonian that forces all the 
chord projections in double phase space to evolve in the same way, independently of 
each centre, Xj(^). Then the evolution of all the semiclassical branches can again be 
united into (15.60 . with the only difference that now $,{t) is obtained from the expansive 
linear equation [14J: 

k = J^ + lt (7.6) 

This is just one of Hamilton's equations for the full double phase space Hamiltonian 
(I6.19p . Since (15.60 coincides with the exact solution of the Lindblad equation for 
quadratic H{x.) and linear L(x), it follows that the present semiclassical theory is exact 
in this limit, even in the presence of dissipation. 
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There is a subtle distinction to be noted in the derivation of the semiclassical 
evolution of a dissipative Markovian system, with respect to the theory in section 5. The 
approximation there could be interpreted as the short time limit of a conceivable periodic 
system, wherein the internal Hamiltonian and the Lindblad operators were alternatively 
turned on and off. The full dissipative approximation is now derived in the same way, 
if we include the new dissipative part in the double Hamiltonian. Even though the 
Trotter ansatz [40j can still be invoked, this mathematical procedure is now devoid of 
its interpretation as a conceivable periodic system, because we cannot physically switch 
off the decohering part of the Lindblad equation without eliminating the dissipative part 
of the Hamitonian: In the present context, it is possible to have decoherence without 
dissipation, but not the other way around. In any case, the extension of the semiclassical 
approximation to dissipative systems is also exact in the quadratic limit. 

The evolution of the semiclassical Wigner function now follows through the 
derivation in section 5: The Fourier transform of (15.21) will be a convolution of Wj^i^, t), 
the Fourier transform of Xj^i^jt)^ with a widening window which coarse-grains over 
the interferences of the Wigner function. Beyond the decoherence time, the Fourier 
transform of the chord function, i.e. the Wigner function becomes smooth and classical- 
like. From then on, the classical motion on this centre space is given by Hamilton's 
(single phase space) equations, with the addition of the purely dissipative term (17.31) : 



(to which diffusion is always added by convolution). The reason is that this is an 
invariant subspace for unitary double phase evolution |25j, a property which is not 
altered by dissipation. Since the decoherence functional has, at this stage, effectively 
cancelled all large chords, the physical interest is concentrated on this plane. In this 
regime the evolution of the decoherenceless factor of the Wigner function in (15. 7p . can 
be pictured as purely classical, iy°j(x, t) = iy°j(x(— t), 0), with x(— t) obtained from 



It might seem strange that the semiclassical solution of the dissipative master 
equation for the chord function becomes exact in the case of linear Lindblad phase 
space functions and a quadratic Hamiltonian. The latter condition is familiar on its 
own, but the exact Lindblad term (I4.10p is not in the same form as the semiclassical 
approximation (14.91) . so that it is not evident that the same limiting behaviour is 
obtained. However, let us, in this case, reinterpret the discrepant term as part of a 
quantized Hamiltonian superoperator for a Schrodinger-like equation corresponding to 
double phase space. Then x and y will be operators corresponding respectively to 
positions and momenta in double phase space, such that x = iH-^. It follows that 
x(— Jy) can then be interpreted as a double wave function, so that 




— 7X 



(7.7) 



(EZD. 



{I' A I") ^ = (r A I") y ■ |- x(-Jy) = -'-{i' A I") y ■ X x(-Jy), 



(7.8) 
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and the action of this Hamiltonian superoperator on a semiclassical branch of the chord 
function, to first order in ft, 



is the same as in (14. 9p . 

8. Semiclassical markovian propagators 

The present semiclassical picture for markovian evolution of the density operator has 
dealt directly with the chord function and its Fourier transform, the Wigner function. 
The problem that must now be addressed is that, whereas the decoherence functional 
quickly quenches the contribution of large chords, the remaining classical small chord 
region lies in the neighborhood of a caustic for both the chord function and the Wigner 
function. In other words, the region of the lagrangian submanifold, A(t), lying close 
to the centre subspace y = in double phase space, projects singularly onto both 
the x-subspace and the y-subspace. Hence, the direct semiclassical theory above is 
only applicable to the initial stages of markovian evolution, before the long chords are 
mostly quenched. Though the further processes of dissipation and diffusion proceed 
continuously through the initial stage and beyond, their direct description requires an 
advanced semiclassical treatment of markovian evolution in the caustic region, beyond 
the scope of the present theory. 

One way out of this problem is to consider alternative lagrangian submanifolds, 
which do not have caustics, at least initially. This is the approach adopted in |25j . 
Instead of the submanifold A(t), which corresponds to pit), but has undesirable caustics, 
we evolve the submanifolds x = constant, or y = constant, corresponding respectively 
to unitary reflection, or translation operators, R^, or T^. The x = constant submanifold 
has no caustic in its y-projection, while the y = constant submanifold is free from 
caustics in its x-projection. Furthermore, a finite time must pass before the evolution 
generated by the double Hamiltonian can bend either of these submanifolds sufficiently 
to produce caustics, until which time both the chord representation of -Rx(^) and the 
centre representation of T^{t) will be represented in the simple semiclassical forms (12.181) 

and (I2.16p . In other words, -Rx(C, t), the chord representation of the reflection operator, 
and T^(x, t), the centre representation of the translation operator, respectively, will 
both have a single semiclassical branch. 

In this way we obtain a semiclassical approximation to the evolution, whether 
unitary or markovian, by inserting the approximate evolved operators in the exact 

ffOn the other hand, it is the essentiaUy nontransversal intersection of the double torus, A = 
£(x^) X £(x+), with the central subspace, along £(x), that allows for a nontrivial Wigner function 
even after decoherence, such as is observed in the quadratic limit. (The single torus, £(x), has the 
same form as the factor tori in the respective spaces, x^.) 




(7.9) 
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relations 



W{^,t) = J ^^x(l)T^(x,t), (8.1) 

or 

xiC t) = j W{^) 2^ R^{t t). (8.2) 

Unlike these mixed propagators, which involve both centres and chords, the semiclassical 
expressions for direct propagators of Wigner functions |17j necessarily involve uniform 
approximations through caustics, even in the simple case of unitary evolution. 

In the present context, it is fl8.2p that should be chosen, because the markovian 
evolution of a reflection operator, initially represented by 



2^^i?.(|)=exp(^-xA^j, (8.3) 

is approximated semiclassically within the chord representation in the same way as 
a single branch of the chord function, x(C) (12.161) . This holds for arbitrary 
Hamiltonians, together with linear Lindblad operators, self-adjoint or not. Thus, the 
approximate evolution of (18. 3p has precisely the same form as (15. 2p . i.e. 

RM.t) = Rlii.t) exp{-l^|*rft'|Z,-|(-t')|'}, (8.4) 

k ^ 

where ^(t) is the trajectory of the initial chord, ^(0) = ^, obtained from the y component 
of the evolving double phase space vector, X = (x, y = J^), and the Lindblad 
coefficients (14.81) are in general complex. Here, the decoherenceless factor is 

RUM = Rl{-Jy,t) = ^|^^i^r''-P (r^^^'^O- ^^-^^ 

The evolution of i?°(^,t) does not coincide with the unitary evolution presented in 
[25] , because of the dissipative term in the Hamiltonian. In short, we merely substitute 
the transported action variables X^, that previously defined the evolving lagrangian 
surface in double phase space, by x in (16.141) and the chord action, crx(^,t) = a'^{y,t) 
is governed by the Hamilton- Jacobi equation (17. 4p . Evidently, the initial action is 
cr^(y,0) = x ■ y, so that the initial semiclassical expression coincides with (18.31) . It 
should be noted that, at the chord origin, R^{$,, t) = for all times, and multiplication 
by the exponential of the decoherence functional does not alter this value. Therefore, 
normalization is preserved, according to (I2.14p . 

The conditions for the derivation of (18. 4p must now be analyzed. Recall that the 
exponential of the decoherence functional was assumed to be a smooth (non-oscillatory) 
function in our previous derivation of the semiclassical chord function (15. 2p . For 
small chords, which are now in focus, this smoothness assumption is falsified as the 
decoherence time is reached, i.e. the time for the volume of the decoherence factor to 
shrink to that of a coherent state. This establishes the duration beyond which this chord 
propagator is valid. Nonetheless, the evolved propagator can be reexpressed in terms of 
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(static) reflection operators, i.e. in the Weyl representation and this whole procedure 
can then be iterated indefinitely. 

It is precisely in this region, where the direct semiclassical approximation for the 
chord function (15. 2p is singular, that the much simpler form, 

2^^°(C,t)=expQx(t)A|), (8.6) 

can be employed [25]. Here, x()f:) is the trajectory, issuing from x(0) = x, that integrates 
(17. 7p in single phase space, because of the restriction to the invariant subspace y = 0. 
This approximation implies that the evolved refiection operator is still represented by 
a vertical subspace in double phase space. This is indeed true for motion generated by 
a quadratic Hamiltonian and even a general closed double Hamiltonian, lHif{X), of the 
form (I6.10p . leads to a submanifold whose tangent space is vertical at the identity 
subspace, throughout the motion, as discussed in [25]. Fortunately, the new open 
dissipative term, 1Hl{X), of the double Hamiltonian (17. ip preserves this feature. In 
essence, the reason for this is that the evolution can be linearized in the neighbourhood 
of a generalized refiection even in the presence of dissipation. 

To evaluate the corresponding approximation for the decoherence functional, we 
expand the double Hamiltonian (I6.19P as 

OH 

i?J(x,y) = — (x)Ay-7x-y, (8.7) 
ax 

which holds to second order in y. Then, Hamilton's equation for the y-motion in the 
neighborhood of the invariant centre subspace is just 

y(t) = [2H2(x(t))J + 7]y, (8.8) 

where the independent centre motion x(t) is determined by the non-hamiltonian classical 
equation (17.71) and H2(x) is the Hessian matrix for iJ(x(t)). The quadratic form 
y • H2(x(t))y can be interpreted as a local Hamiltonian for the motion transverse to the 
centre subspace, though it receives an extra boost from the dissipation coefficient, 7. 
Given an initial point X = (x, y = J^), whose orbit is assumed to remain close to the 
identity subspace within the time t, we thus obtain its chord evolution as 

^,(x,0 = G,(x)|, (8.9) 

where the classical propagation matrix is 

^ rt nt 1 

G,(x) := lun^ n [77i2H2(x(-))J + 7]], (8.10) 

which, for small times is approximately 

Gt(x) ^ exp [ /* dt' (2H2(x(t')) J + 7)1 • (8.11) 

L JO J 

The only difference between this approximation and the exact evolution for the quadratic 
case in [13] is the x-dependence for H2 and, hence, the need for an integral in the 
definition of G^. Furthermore, the present definition incorporates the dissipation 
coefficient, 7, from the full double Hamiltonian (16.191) . 



Semiclassical Evolution of Dissipative Markovian Systems 



28 



We are dealing here with trajectories that remain very close to the invariant 
centre subspace, where D{x(t), y = 0} = 0. In other words, within the short chord 
approximation, the pair of trajectories, which define the decoherence functional, become 
indistinguishable from a single classical trajectory, so that 

fl{xW.y-0}=^^. (8.12) 

a simple quadratic form in the variables transverse to this centre subspace, specified by 
the evolving matrix, 

Mi(x) = fdt' Gnx)Uf G^x). (8.13) 

Here (.)"^ denotes the transpose of a matrix, or a vector. This quenching exponent has 
exactly the same form as the exact one for quadratic Hamiltonians [H] along the centre 
subspace. The only difference lies in the dependence of the classical chord propagator 
Gi(x) on the local quadratic approximations of the Hamiltonian along the trajectory 
x(t), which can be taken either in the forward direction, starting at x, or backwards 
from x(t). In conclusion, the small chord approximation for the mixed propagator is 

2^ Utt) = exp (jx(t) A ^) exp ( - ^^^^). (8-14) 
This can now be inserted into (18.21) . to obtain 

Here, we cannot use the semiclassical approximation (12.181) for the initial Wigner 
function, Vr(x), in the above integral, because it is singular in the small chord region, 
but improved uniform approximations [32l HH] are valid in this range. A correct initial 
normalization of the Wigner funtion guarantees the normalization of the chord function 
for any subsequent time, according to (12.141) . 

It is now worthwhile to review the full construction required to follow the markovian 
evolution of a given initial pure state. By inserting its Wigner function into (18.21) along 
with the semiclassically evolved propagator, specified by (18. 4p and (18.51) . we obtain 
the evolution of the chord function for a time interval of the order of the decoherence 
time. The amplitudes of all long chords will be strongly quenched in this evolved 
chord function, so that its Fourier transform, the evolved Wigner function, becomes 
smooth and positive. This can be further evolved by inserting it into (18.21) again. But 
now there will be no large chord contribution, even right from the start, because at 
t = 0, equation (18.21) reduces to the inverse of the Fourier transform that we have just 
made. Therefore, the much simpler evolution given by the small chord approximation 
(I8.15P is now adequate. The markovian evolution never dequenches the large chord 
amplitudes, so that all further iterations of this procedure may safely rely on the small 
chord approximation. 

Let us now describe the evolutions that can be calculated within the small chord 
approximation entirely within the Weyl representation. In the limit of small propagation 
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time, M( 0, the Fourier transform of the evolving chord function is exactly 
W{x',t) = iy(x(— 1),0), the classical propagation of the Wigner function. This is 
the limit where the Wigner- Wigner propagator becomes merely 5(x' — x(t)). Unless 
the Hamiltonian is quadratic, it is only valid because the long chords that have already 
been quenched out of the propagator. For longer times, it is also possible to perform 
the Fourier transformation of flS.lSp exactly, so as to obtain the Wigner function as an 
evolving convolution with a gaussian that broadens from an initial Dirac (5-function. 
Following [Hj, this is given by 

Wt^'t\-r "^^ .xnl' (x-x-(-t)).M;(x)(x-x-(-f)) , 



in which x'(— t) is the backward trajectory of the evaluation point, x'. It is important to 
emphasize here, once again, that this is a trajectory not of the single Hamiltonian H{x.), 
but of equation (17.71) . which adds the dissipative term, — 7x, to Hamilton's equation for 
the single Hamiltonian -ff(x). The matrix 

M;(x) := -JMt(x)-ij (8.17) 

characterizes the Fourier transform of the exponential of the decoherence functional in 
f lS.lSp . This is now a broadening gaussian window, which coarse-grains the classical 
evolution of the initial Wigner function. The same coarse-graining, which accounted for 
the initial loss of quantum coherence of an initial pure state, can now be interpreted 
as resulting from a classical Langevin equation for Brownian motion acting on a purely 
classical probability distribution in phase space. This description is again the exact 
result in the quadratic case [H] and it has been derived independently many times 
for the dampened harmonic oscillator ^151 [T6l [T7] . In the purely classical context, this 
picture for Brownian motion can already be found in the review |49] . 

It is tempting to extrapolate the approximate small chord evolution (18.151) beyond 
the decoherence time. After all, the effect of this process is precisely to eliminate the 
contribution of long chords, so that the passage to (I8.16P becomes more valid. Indeed, 
there would be no obvious contradiction if the validity of (I8.15P were to extend beyond 
the range allowed by our derivation, as is true in the quadratic case. If we do insert 
(I8.14P into (14. 4p . the integral to be evaluated becomes 

{LR^P - h^LR^ - ^R^Pl){C) = 

{l(x"+ ^W- |^)*-^[L(x"+ |)* + L{^"- ^W- y)*]}- (8.18) 

Recalling that the Lindblad functions i^(x) are assumed linear, we can now substitute 
x" by iHJd/d^' within the brackets {} in (I8.18P and then integrate over this variable to 
obtain 6{$,' — ^"). Thus, in this case, the integral of (I8.18P becomes 

{LR^D - h^LR^ - ^R^DL)i^') = 
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A r')[x(t) A r - 1' • M,(x)^'] - ■ a' + H" ■ O']) Mi\ t) (8.19) 

if we neglect semiclassically small terms arising from the derivatives of the Lindblad 
functions. 

This is almost in the same form as (14.91) . except for the term A I") ■ Mt(x)^', 
which is not semiclassically small. Nonetheless, this term does disappear if the 
dissipation coefficient 7 = (Z' A I") = 0, so allowing us to carry through an analogous 
derivation to that of the semiclassical chord function in section 5. The result is that 
in this case the small chord approximation (I8.15P is valid and the closed formula for 
long term evolution of the Wigner function (18.161) can be legitimately extrapolated 
far beyond the initial period for decoherence. In the general dissipative case with 
nonquadratic Hamiltonians, the qualitative picture of the Markovian evolution, resulting 
from repeated iterations of (18.161) is the same, but a single closed formula is not yet 
available. 

9. Discussions 

The point of view of this paper is that semiclassical Wigner functions or chord 
functions stand to double phase space as do semiclassical position and momentum wave 
functions to single phase space. In strict analogy to the more familiar theory, each 
of these conjugate representations is defined in its own subspace and contains complete 
information concerning the quantum state, be it pure or mixed. However, mixed systems 
demand a density operator description, rather than as states in Hilbert space, which 
can be provided by the Wigner function or the chord function. For strictly unitary 
evolution, it is still possible to restrict consideration to the single phase spaces on which 
these functions are originally defined, but not for general quantum markovian processes. 

No attempt has here been made to expand the framework of quantum mechanics 
itself. The Wigner function and the chord function are particular choices of 
representation for the evolving density operator, each in their own phase space. It 
is only the identification of the various terms in these representations of the Lindblad 
equation with those of standard semiclassical wave function evolution, that leads to a 
generalized WKB-like solution for quantum markovian motion. 

The full double phase space Hamiltonian would have an unfamiliar form if it were 
to be considered as the generator of motion for a mechanical system in an ordinary 
multidimensional phase space, but, once these peculiarities are understood, there results 
a qualitative picture for generalized semiclassical evolution that is pleasingly intuitive: 
The decoherence functional quickly quenches the contribution of all large chords, just as 
in the exact quadratic case [H]. Hence, after a short decoherence time, we may restrict 
the analysis to the neighbourhood of the centre subspace y = in double phase space. 
The fact that this is the subspace, where the Wigner function is supported, indicates 
that the conjugate pair of the Weyl and the chord representations constitute a privileged 
frame for the study of markovian evolution for quantum systems. 
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The various semiclassical representations of the density operator are derived from 
a single Lagrangian submanifold in double phase space, with its hamiltonian evolution. 
It so happens that the product nature of this manifold implies that its intersection 
with the identity subspace coincides with a caustic for both the Wigner function and 
the chord function. We have shown in section 8 that a simple adaptation of the 
present theory, obtained by choosing a set of more favourable lagrangian submanifolds, 
furnishes an optimum semiclassical propagator. This allows for appropriate small chord 
approximations, leading back to a closed formula for the evolution of the Wigner function 
f l8.16p . Thus, we again establish contact with the exact markovian theory for quadratic 
Hamiltonians. 

Notwithstanding that the evolution of our propagator is derived in the semiclassical 
approximation, it will transport any kind of pure or mixed Wigner function. These may 
be initial semiclassical states, but also (squeezed) coherent states, Schrddinger cat states, 
or whatever. Though the analysis is harder in the intermediate stage, between the initial 
loss of quantum coherence and asymptotic classical motion, it can be conjectured that 
the Wigner function becomes positive everywhere, even if it is doubtful that the time 
for this will be independent of the initial state, as in the quadratic case [SI [Hj. 

The fact that the present semiclassical theory is exact in the case of a quadratic 
internal Hamiltonian, even in the presence of linear Lindblad operators, can be 
considered as an indication that it provides a useful generalization of this simple case. It 
will certainly be necessary to make detailed comparisons of the approximate semiclassical 
evolution to the direct integration of the exact equation in the case of nonquadratic 
Hamiltonians. However, the integration of the multidimensional partial differential 
master equation is a considerable enterprise, specially for the highly oscillatory Wigner 
functions described by semiclassical theory. An alternative is to resort to quantum monte 
carlo methods, such as in [3], but then the comparison is merely between alternative 
approximat ions . 

Our analysis has dealt only with Lindlad operators that are linear functions of 
position and momentum operators. This may be justified by picturing these operators 
as the quantum variables that are responsible for the coupling to the environment, as 
in the derivations of our standard example ( 11. 2p . Weak coupling, as assumed for a 
markovian theory, often implies that a linear approximation is valid, but this need not 
be so. Perhaps it has been the very difficulty of obtaining a fully reliable general picture 
of the evolution of the density operator of an open system that has so far hampered the 
study of systems described by nonlinear Lindblad operators. Further generalization of 
the present theory to include this possibility is the subject of ongoing work. 

It is early to predict whether the insight provided within a full semiclassical theory 
of non-unitary evolution will reveal features which are beyond our present intuition. This 
can only be decided by analyzing examples of increasing complexity within the present 
framework, or by answering harder mathematical questions of the formalism itself (e.g. 
as in [50]), or perhaps both. For the present, we have gained in understanding how 
dissipation and diffusion are emmerging properties of the single continuous markovian 
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evolution in double phase space, which, in its initial stages, is best described as 
decoherence. It is only when the contribution of long chords is quenched onto the 
neighbourhood of the centre subspace that the motion can meaningfully be described 
as dissipative. 
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